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Abstract 



We study the contribution of the gluino-squark mediated flavor changing process for the 
CP violation in 6 — )• s and 6 — t- d transitions facing on recent experimental data. The 
mass insertion parameters of squarks are constrained by the branching ratios of 6 —t- 57 
and 6 — 7- (i7 decays. In addition, the time dependent CP asymmetries of — )■ (j)Ks and 
—7- rj'K^ decays severely restrict the allowed region of the mass insertion parameter 
for the b ^ s transition. By using these constraints with squark and gluino masses of 
1.5 TeV, we predict the CP asymmetries of Bs — s- (jxp, Bg — )■ rj'cj), and B^ — t- K^K^ decays, 
as well as the CP asymmetries in 6 — t- S7 and b ^ d'j decays. The CP violation in the 
Bs — )• <j)(j) decay is expected to be large owing to the squark flavor mixing, which will be 
tested at LHCb soon. 
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1 Introduction 



The LHC experiments are now going on to discover new physics, for which the supersymmetry 
(SUSY) is one of the most attractive candidates. The SUSY signals have not yet observed 
although the Higgs-like events are almost confirmed [1]. The lower bounds of the superparticle 
masses increase gradually. The squark and the gluino masses are expected to be larger than 
1 TeV [2]. 

On the other hand, the LHCb collaboration has reported new data of the CP violation of 
B mesons and the branching ratios of rare B decays The new physics is also expected 

to be indirectly found in the B meson decays. For a long years the CP violation in the K 
and B^ mesons has been successfully understood within the framework of the standard model 
(SM), so called Kobayashi-Maskawa (KM) model [6], where the source of the CP violation 
is the KM phase in the quark sector with three families. While, there are new sources of 
the CP violation if the SM is extended to the SUSY models. The soft squark mass matrices 
contain the CP-violating phases, which contribute to the flavor changing neutral current 
(FCNC) with the CP violation. Therefore, we expect the effect of the SUSY contribution 
in the CP-violating phenomena. However, the clear deviation from the prediction of the SM 
has not been observed yet in the LHCb experiment [31 IH [5]. 

In our previous works [3 [S] , we studied the SUSY contribution, which comes from the 
gluino-squark mediated flavor changing process We used the only experimental data 

of the 6 — S7 decay to constrain the mass insertion (MI) parameters of squarks. And then, 
we predicted the CP violations of a few 6 — t- s transition processes. In the present paper, 
we give the systematic studies for the effect of the gluino-squark mediated flavor changing 
process in the CP violation of the — >■ s and b ^ d transitions. In order to obtain more 
precise numerical results, we take account of the QCD corrections for the SUSY contribution. 
Moreover, in order to constrain the MI parameters, we also input the recent experimental 
data, the time dependent CP asymmetries of B^ non-leptonic decays, in addition to the 
experimental data of the b S'j decay. 

The LHCb collaboration reported the time dependent CP asymmetry Sj/^^ in the non- 
leptonic Bs decay, which gives a constraint of the SUSY contribution on the 6 — )■ s transition. 
The CP asymmetry of Bg — (fxp is expected to be observed in the near future at LHCb [S] . 
If there is the contribution of the squark flavor mixing in the FCNC, we expect to observe 
the sizeable time dependent CP asymmetry in this process, in which the SM prediction is 
very small. 

The typical process of the b s transition is the b ^ decay, in which the experimental 
data of the branching ratio, the direct CP violation, and the time dependent CP asymmetry 
Sk*'-! have been reported. The SUSY contribution is also constrained by the data of the time 
dependent CP asymmetries in B^ — j- (pKs and 5° — rj'K'^ decays [221 123] . 

On the other hand, the b ^ d transition also becomes available to investigate the SUSY 
contribution quantitatively taking account of the recent experimental branching ratio of 
the b ^ d'-f decay [211 125] • In this transition, the time dependent CP asymmetry of the 
B^ —7- K^K^ decay is an attractive one to search for the SUSY effect because the penguin 
amplitude dominates this process. We also predict the time dependent CP asymmetry of 
5° ^ p7, Sp,. 



1 



The dominant contribution of the SUSY is the gluino-squark mediated flavor changing 
process for the B meson decays discussing in this work. We present the constraint for the 
MI parameters {6^^)23 and (5^'^) 13 by putting the experimental data, where squarks and 
the gluino masses are at the TeV scale. By using these MI parameters, we predict the CP 
violation of the B meson decays, in which the interesting one is the Bg — )■ (jxp decay. The CP 
violation in this decay will be measured at LHCb in the near future. 

In section 2, we present the formulation of the gluino-squark contribution on the CP 
violation of B mesons in our framework. In section 3, we discuss the 6 — )■ s transition, and 
present numerical predictions for the direct CP violation and time dependent CP asymmetries 
in 5° — >■ K*'j, Bg — )■ 00, and Bg — >■ rj'cj) decays. In section 4, we discuss the b ^ d 
transition, and present numerical predictions for the direct CP violation and time dependent 
CP asymmetries in B^ — p7 and B^ — )■ K^K^ decays. Section 5 is devoted to the summary. 



2 Squark flavor mixing in CP violation of B mesons 



Let us present the framework of the calculations for the contribution of the squark flavor 
mixing, which is the coupling among down-type quarks, down-type squarks, and the gluino. 
The effective Hamiltonian for the Ai? = 1 process is given as 



H, 



71" 



q'=u,c i=l,2 

where q = s,d. The local operators are given as 



i=3-6,77,8G 



O3 = iqal^PM Y^iQp^PLQp), O4 = {qal^Pibp) Y,{Qp1^PlQo), 

Q Q 

O5 = {qalf^PM Y.^Qp1^PrQp). Oe = {q^-f^PLbp) Y,{Q pl^ PrQ , 



O 



77 



jmbq^a'^'' Pjib 

a Pi 



o. 



8G 



;m,q^a^''PRT\bf,G^ 



(2) 



where Pr = (1 + 75)/2, Pl = {1 — 75)/2, and a, (3 are color indices, and Q is taken to be 
M, (i, s, c quarks. Here, Cj's and Cj's are the Wilson coefficients, and Oj's are the operators by 
replacing L{R) with R{L) in Oj. In this paper, Cj includes both SM contribution and gluino 
one, such as Ci = Cf^ + Cf, where Cf^'s are given in Ref. [26] . 

In order to estimate the SUSY contribution for Cf, we take the most popular ansatz, a 
degenerate SUSY breaking mass spectrum for the flavor structure of squarks. In the super- 
CKM basis, we can parametrize the soft scalar masses squared of the down-type squarks. 
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M? , M? , M? , and M? as follows: 

a-LL fl/ifi OLfl QflL 



d.LL 



Mi 



Mi 



771: 



mis 



11 



13 



1 + (^i'')33. 

(5^)13 
22 )23 
(5^)23 1 + ('5^)33. 



112 

1 + (^i^)22 
('5i")23 

1 + (5^^ 



(M; 



2 )t 



(^r)i2 (^i^)i3 



21 (<^d^)22 {Sd^)23 
31 (<^d )32 (<^d )33; 



(3) 



where is the average squark mass, and and are called as 

the mass insertion (MI) parameters. 

The Wilson coefficients of the gluino contribution Cf are given as follows 



Cf(m,) 



9 
1 



10 ^ , , 



-Bo(x) 



18 ^ ^ 
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:P2(x) 



-P^ix) 



— Pl(x) - -P2(X) 

18 ^ ^ 2 ^ ^ 
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Plix) 



:P2(x) 



6GFVt,V,lml 



fc3 



-M3(x) -/itan/3-|-M,(x) + {6ji'^)k3—-Mi{x 



m? 3 



nib 3 




-Msix) + 3M^{x) 



/itan/3— I ( -M„(x) + SMJx 
mi V 3 




;fc3- 



1 



-M,{x) + 3M2{x] 



(4) 



where /c = 2, 1 correspond to b ^ q [q = s,d) transitions, respectively. Here the double mass 
insertion is included in Cj^{mg) and CgQ{mg). The Wilson coefficients Cf{mgys are obtained 
by replacing L{R) with R{L) in G^{'mg)^s. The loop functions in Eq.(|l]) are presented in our 
previous work [7]. In our calculations, Cjj and Csg give dominant contributions to the CP 
violations in 6 — s and b d transitions. The effective Wilson coefficients of C-j^{m}j) and 
CsG{f^b) are given at the leading order of QCD as follows [26] : 



Csci^b) = vCiG{m-g), 



m 



'g>i 



(5) 



where 

16 \ ^ \ — \ — 



Let us discuss the time dependent CP asymmetries of and Bg decaying into the final 
state /, which are defined as [28] 



2IniA/ _ 1 - |A/ 

where 
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Here M^2(5' = c?) are the dispersive parts of the Bq-Bq mixing, in which the quark-squark- 
gluino interaction contributes in addition to the SM one. The MI parameters {5^^)kz and 
{S^^)k3 {k = 2, 1) are constrained by CP violations in the AB = 2 transition as discussed in 
our previous works [71 IH]. On the other hand, (5rf^)fc3 and {S^'")k3 {k = 2, 1) are constrained 
in the AB = 1 transition. 

In the 5° —7- J /ipKs and Bg — )■ J /'>p(i) decays, we write Xj/^Ks and Xj/^^ in terms of phase 
factors, respectively: 

Aj/^^, = -e-^'^^ Xj/^^ = e-'^\ (9) 

The recent experimental data of these phases are [231 EH] 

sin <pd = 0.679 ± 0.020 , 0, = -0.002 ± 0.083 ± 0.027 . (10) 

We expect the SUSY contribution to be included in these observed values. 

Since the B^ — )■ J/ipKs process occurs at the tree level in the SM, the CP asymmetry 
mainly originates from Mf2- Although the B^ — )■ (pKs and B^ — )■ rj'K^ decays are penguin 
dominant ones, their CP asymmetries also come from in the SM. Then, the CP asym- 
metries of B^ — > J/ipKs, B^ —7- (f)Ks, and B^ — )■ rj'K^ decays are expected to be the same 
magnitude. 

On the other hand, if the squark fiavor mixing contributes to the decay at the one- 
loop level, its magnitude could be comparable to the SM penguin one in B^ — )■ (pKs and 
B^ — )■ rj'K^ decays, but the squark fiavor mixing contribution is tiny in the B^ — )■ J/ipKs 
decay because this process is at the tree level in the SM. Therefore, there is a possibility to find 
the SUSY contribution by observing the different CP asymmetries among those processes [30] . 

The time dependent CP asymmetry Sj/^Ks has been precisely measured. On the other 
hand, PDG [22] and Heavy Flavor Averaging Group (HFAG) [23] presented considerably 
different values for S^Ks while almost same one for S^'ro- Each of the observed ones in 
HFAG is consistent with the SM prediction. In order to get conservative constraints, we take 
the data of these time dependent CP asymmetries in HFAG [23], which are 

Sj/^Ks = 0.679 ± 0.020 , S^Ks = 0.74l°;JJ , 5^/^o = 0.59 ± 0.07 . (11) 

These values may be regarded to be same within the experimental error bar. Thus, the 
experimental values are consistent with the prediction of the SM. In other words, these data 
severely constrain the MI parameters (5^^)23 and {5^^)23 in our following analyses. 
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3 The 6 — > s transition 



At first we discuss the contributions of tlie squark flavor mixing for the 6 — )■ s transition, 
which are given in terms of the MI parameters (5^'^) 23, (5j^^)23, (5^'^) 23, and (5^^)23. These 
MI parameters are constrained by the experimental data of B meson decays. 

Let us show the formulation of the 6 — s transition. The CP asymmetries Sf of Eq. ([7j) 
in the b — sss transition are one of the most important processes when we investigate the 
new physics. The CP asymmetries Sf for — )■ (pKs and B^ rj'K^ are given in terms of 
A/ in Eq. 

J2 (Cf^(0.) + Cf(0.)+Cf(0.)) 

\ i0 i=3-6,77,8G , - 

j=3-6,77,8G 

where (O,) is the abbreviation of (/ 10,15°). It is noticed {(l)Ks\Oi\B°) = {(f)Ks\di\B°) 
and {ri'K^\Oi\B^) = —{ri'K^\Oi\B^), because these final states have different parities [30| [3T]. 
Since the dominant term comes from the gluon penguin C|g, the decay amplitudes of / = (pKs 
and / = rj'K^ are given as follows: 

1(5° ^ (j)Ks) oc CsGinib) + CsGinib), 

A{B' ^ V^O) oc Cscim,) - Cscim,). (13) 

Since Csci'm'b) is suppressed compared to Csci^b) in the SM, the magnitudes of the time 
dependent CP asymmetries Sf (/ = J/ip4', (pKs, rj'K^) are almost same in the SM prediction. 
However, the squark flavor mixing gives the unsuppressed Csci'^b), then, the CP asymmetries 
in those decays are expected to be deviated among them. Therefore, those experimental data 
give us the tight constraint for Csci^b) and CsciiTT-b)- 
We have also A/ for Bs 00 and Bg — )■ (pf]' as follow: 

Xu H = e-'^' -^Z^-^ , (14) 

''''' Y: cr*m+crm+crm 

i=3-6,77,8G 

with (00|Oj|i?s) = — (00|Oj|-Bs) and (0?7'|Oi|i?s) = (0?7'|Oi|i?s). The decay amplitudes of 
f = (j)(f) and / = (prj' are given as follows: 

A{Bs 00) oc Cscimb) - Cscimb), 

A{B, ^ 0r/') cx Cscimb) + Cscimb). (15) 

Since CsciOsc) and Csg{Osg) dominate these amplitudes, our numerical results are insensi- 
tive to the hadronic matrix elements. In order to obtain precise results, we also take account of 
the small contributions from other Wilson coefficients Ci [i = 3, 4, 5, 6) and Ci [i = 3, 4, 5, 6) 
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in our calculations. We estimate each hadronic matrix element by using the factorization 
relations in Ref. [32]: 



{Osg) 




N, 



(O4) + (Oe) - — ((O3) + (O5)) 



(16) 



where (g^) = 6.3 GeV^ and A'^c = 3 is the number of colors. One may worry about the reli- 
ability of these naive factorization relations. However, this approximation has been justified 
numerically in the relevant b —> s transition as seen in the calculation of PQCD 

Let us discuss the contribution of the MI parameters to C^q in the Eq. (jlj). Since the 
loop functions are of same order and rriq ~ rrig, the ratio of the LL component and the 
LR one is {6^^)23 x yutan/^/m^ to {6^^)23 x rriq/mb. If C(/x tan/3) ^ 0{mq) and > 1 
TeV, the LR component may contribute significantly to C|g due to the enhancement factor 
mg/rrib = C(IO^). For example, in the case of (5^'^)23 = 10"^ and {8^^)23 = 10"^, the 
LR component dominates C^q, while it is minor in Mf^ [HIH]. Actually, the magnitude of 
{8d^)23 is at most 10~^, which was estimated in our previous works [3 [8]. In our following 
calculations, we take |(5^'^)23| ^ 10~^. 

We can also constrain the SUSY contribution from the 6 — 57 decay. Here we discuss 
three observable values, those are the branching ratio BR(6 — 57), the direct CP asymmetry 
A^"^*^, and the time dependent CP asymmetry of — )■ A'*7, Sk'-j- The branching ratio 
BR(6 — g7)(g = s, d) is a typical process to investigate the new physics. It is given as 



6a 



BR(6 ^ g7) 
BR(6 ^ cez7,) iV^bP nf{z) 



(|C7^(mfe)|2 + |C'7^(mb)p 



(17) 



where 



\2z^\^z 



m 



c,pole 



m 



b,pole 



Here Cj^{mij) and Cj^{m}j) include both contributions from the SM and the gluino-squark 
mediated flavor changing process at the mb scale. As seen in Eq. MI parameters (5^^)fc3 
dominate both Cj^ and Cf^. Therefore, we should discuss the contribution from {5^^)k3 in 
our numerical calculations. 

We can also estimate the direct CP violation A^^^^ b ^ q'j decay (g = s,d), which 

is given as [35] 



_ r{B ^ x,7) - r(fi ^ x,-7) 
r(5 ^ x,7) + r(5 ^ x,-7) 



£;t,>(i-<5)£;" 



^Im[C2Cy-|[^(.) + 6(^,5)]lm 



1^ ^ ^^^^^ 



l^2'^8G 



(19) 
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where v{z) and b{z,6) are explicitly given in Ref . [35] . and Ci, Ci {i = 77, 8G) include both 
the SM and SUSY contributions at the mf, scale. 

The time dependent CP asymmetry Sk'-j in the — )■ K*'-/ decay is also important 
measure of the CP violation: 

\C,,{m,)/C-r,{m,W + l 

This CP violation comes from the interference between Cj^{mi,) and Cj.~f{m}j) [271 ES]- In 
the SM, C^^{mh) /Cj}:l^{mij) oc ms/mi, for this process. Therefore, Sk*'^ is suppressed [36] . 
However, Sk*'^ could be enhanced owing to the squark flavor mixing. 

Our setup in our calculations are shown as follows. We take tan /3 to be 1 TeV, and set 
|(^i''^)23| — |(<^d^^)23| ^ 10~^ following from our previous works [71|8]. Then, the contribution 
of these MI parameters to C^^ and C^q are minor. On the other hand, (5^'^)23 and (5j^^)23 
are severely constrained by magnitudes of Cj^ and C^g- In addition, we suppose \{5^^)22,\ = 
l('^d^^)23|- Then, we can parametrize the MI parameters as follows: 

(5.^^)23 = |(5i^)23|e^^^- , (5^)23 = |(5."^)23|e-^- . (21) 

Now we show numerical analysis in our setup. In our following numerical calculations, 
we fix the squark mass and the gluino mass as 

mq = 1.5 TeV, = 1.5 TeV, (22) 

which are consistent with recent lower bound of these masses at LHC [2]. 

In our analysis, the present experimental data of BR(6 — 57), Sj/^Ks^ S^Ksi and S'^'a'o 
give tight constraints for MI parameters. Here we put the experimental data [22] 

BR(6 S7)(exp) = (3.53 ± 0.24) x 10-^ (23) 

on the other hand, the SM has predicted [37] 

BR(6 ^ S7)(SM) = (3.15 ± 0.23) x 10"^ (24) 

Therefore, there is a room for the contribution of the gluino-squark mediated flavor changing 
process^ For Sj/^Ks-, S<t>Ks^ and Sn'RO, we put the data in Eq. lfTTl) . In the SM, these 
magnitudes of Sj agree with among them. 

At first, in Fig. [1] (a), we show the allowed region in the plane of the absolute value 
I ((5^^)23 1 and the phase , where the only experimental constraint of BR(6 — 57) is put. 
The magnitude of the MI parameter (5^*^)23 is allowed as | (5^^)23! ^ 9 ^ 10~^. We note 
that the SUSY contribution to BR(6 — )■ 57) becomes superior compared with the SM in 
the region of |((5^'^)23| ^ 4 x 10~^. It is also noted that the {6^^)23 is almost real around 
the upper bound 9 x 10~^, that is at 26*23^ = or 27r. The experimental constraints of 
Sj/^Ksy S(i,Ks, and S^^'ko give the severe cut as seen in Fig. [1] (b), where these experimental 
data are put in addition to BR(6 — > 57). In this figure, any value of the phase is allowed 

^In our analysis, we do not take account of the contribution of the charged Higgs and chargino in 6 — ^ 57. 
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in 1(5^^)23! ^ 5 X 10 ^. On the other hand, the larger region of |(5^'^)23| is allowed until 
2 X 10-2 around the specific 9^^, n/A and Stt/aE The obtained bound 1(5^^)231 < 2 x 10"^ 
depends on the gluino and the squark masses. If they increase, the upper bound is rescaled 
approximately as |(5^'^)23| x 'mq/{1.5 TeV). 

By using this allowed region of {6^^)23, we predict A^(^^^ , Sk'^, S^j,^^ and 5*0^/. In Fig. [21 we 
show the predicted direct CP asymmetry A!^^^"' versus |((5|'^)23|- Here the value at |(5|'^)23| = 
is the SM one, Acp''^(SM) ~ 4 x lO'^ [35]. We predict -3 x lO'^ < A^p'^ < 3 x 10"^ owing 
to the squark flavor mixing. Recent experimental data is still consistent with our prediction 
due to the large error as seen in ^^'^^'''(exp) = —0.008 ±0.029 |22]- The precise data will give 
us an additional constraint of the MI parameters in the future. 

In Fig. [31 we show the predicted CP asymmetry, Sk*'^ The predicted value in the SM is 
S'/^.^(SM) ~ (2ms/mb) sin^^ ~ 4 x 10"^ [36], while the experimental result is S'x*7(exp) = 
—0.15 ± 0.22 [22]. Our prediction is —0.4 < 5";^.^ < 0.2, which is still consistent with the 
experimental data. We also expect the precise data in the near future to test our prediction. 

Although the experimental data of the time dependent CP asymmetries S^Ks ^^nd S^t^^ 
are taken as the input in our analysis, these calculated values do not always cover all ex- 
perimental allowed regions due to the constraint from BR(6 — 37). Those allowed regions 
are shown in Fig. [H The SM prediction is Sj/^i^g (SM) = S'^i^g(SM) = 5'^/ii'o(SM), while the 
present data of these time dependent CP asymmetries are given in Eq. ( ITT]) . The region of 
the right-down corner in the figure is excluded. It is testable in the future experiments. 

In Fig. [5l we predict the time dependent CP asymmetries S^^ and S^ri'- These CP 
asymmetries must be equal to Sj/^^ in the SM. We use the experimental result of S'j/^,^ for 
the phase 0s, which is given in Eq. (fTOj) . in our calculations. We denote the small green line 
as the SM value S'j/^^(SM) = — 0.0363tQ;oQ}5 [M] in the figure. In conclusion, we predict 
—0.2 < S^^ < 0.4 and —0.5 < S^ri' ^ 0.4, respectively. Since the phase 0^ has still large 
experimental error bar, our prediction will be improved if the precise experimental data of 
Sji-^cj) will be given in the near future at LHCb. Since the time dependent CP asymmetry 
will be measured at LHCb, our prediction will be tested soon. 



^Thcrc still remains a very small allowed region around |((S^^)23| = 9 x 10~^, where {8^^)23 is almost 
real, since this region cannot be excluded by the time dependent CP asymmetries. In our work, we omit this 
region hereafter. This region is uninteresting because the SUSY contribution is much larger than the SM one 
in & — S7. 
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(a) 



(b) 




I(<5r)23|xl0^ 



^ 7r/2 




0.5 1.0 1.5 

I((5f)23|xl02 



2.5 



Figure 1: The predicted region of (5^ )23- In both figures (a) and (b), the horizontal and 
vertical axes denote the absolute value and the phase of {6^^)23, respectively. In the figure 
(a), the only experimental constraint of BR(6 — )■ 57) is taken account. In the figure (b), the 
experimental constraints of BR(6 — )■ 57), S^^Ksy and Sr,'KO are taken account. 




-0.06- 



0.0 0.5 1.0 1.5 2.0 2.5 

I(<5!/'')23|X102 

Figure 2: The predicted direct CP asym- 
metry A'^(^^'^ of 6 — s7 versus | (5^^)23!. 
The red solid and two red dotted lines de- 
note the best fit value, upper and lower 
bounds of the experimental data with 90% 
C.L., respectively. 
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Figure 3: The predicted CP asymmetry, 
Sk*^ of 5° i^*7 versus | (5^^)23!, where 
the red solid and two red dotted lines de- 
note the best fit value, upper and lower 
bounds of the experimental data with 90% 
C.L., respectively. 



4 The h ^ d transition 

In this section, we discuss the h ^ d transition as the same way in the b ^ s one. The 
SUSY contribution is given in terms of the MI parameters (^^■^)i3, {Sd^)i3, {^d^)i3, and 
i^d^)i3- The typical b —> d transition is the b ^ d'y decay. The experimental data of its 
branching ratio gives the constraint for these MI parameters. By using these MI parameters, 
we calculate SUSY contributions to the direct CP violation of the 6 — > ^7 decay and the 
time dependent CP asymmetry in the — )■ p7 decay. We also predict the time dependent 
CP asymmetry of the — )■ K^K^ decay. 
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0.4 ^ 1 

0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Figure 4: The allowed region of the time 
dependent CP asymmetries on the S^Ks^ 
Sn'K" plane. The SM prediction Sj/^Ks = 
S<i>Ks = Sri'RO is plotted by the green slant 
line. The experimental data with error 
bar is plotted by the red solid hnes at 90% 
C.L.. 
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-0.6 ^ ^ 

-0.4 -0.2 0.0 0.2 0.4 0.6 

S(f,(i, 

Figure 5: The predicted time dependent 
CP asymmetries on the S^^-S^r)' plane. 
The small green line denotes the SM pre- 
diction from the experimental data of 



In order to constrain the MI parameters, we input the experimental data of the branching 
ratio of 6^ ^7 [211 [25], 

BR(6 ^ ci7)(exp) = (1.41 ± 0.57) x 10-^ (25) 

on the other hand, the SM has predicted [25] 

BR(6 ^ rf7)(SM) = {IMtofi) X 10^^ (26) 

Next we present the formulations of the time dependent CP asymmetries and direct CP 
violation including SUSY contributions. The branching ratio and direct CP violation in 
the b ^ dj decay are given in Eqs. ( |T71) and ( fT9|) . respectively. The time dependent CP 
asymmetry Sp.y in the — )■ p7 decay is an important observable to search for the new 
physics and given as 

Since Cf^{'mh)/Cf^{mb) oc md/rrib in the SM, Sp^ may be expected to be quite sup- 
pressed [3H]- However, Sp^^ could be also enhanced owing to the gluino-squark mediated 
flavor changing process. 

The time dependent CP asymmetries Sj^o^o and Cj^o^o in the — j- K^K^ decay are 
also interesting ones to search for the new physics since there is no tree process of the SM in 
the 5° K^K^ decay [Ml SO]. These CP asymmetries are given in Eq. (JTj) as 

2ImA^o^o 1 - \\koko? /r,oN 

^K^RO — 1 , I X ; '^K^K^^ — -1,1, ; l^Oj 
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where 



\koko 



-P , 
P 



1 
P 



(29) 



The amphtude A(B^ — )• R^R^) is given in Ref. [39], I in which the QCD corrections are 
important for the hadronic matrix elements [H], as 



4Gf 



q=u,c 



Here X is the factorized matrix element (See Ref. [39].) as 

X = -ifKFQ{m\){m\ - m\) 



(30) 



(31) 



where fx and Fo(m^) denote the decay coupling constant of the R meson and the form 
factor, respectively, and = 2m\/{{mi, — ms){ms + rrid)) denotes the chiral enhancement 
factor. The coefficients a?'s are given as 



al{mb) = (C4 - C4) + 



(C3 



C3) ^ as{mb) Cf 



al{mb) = (Ce - C,) + 



(C5 



(C4 - d) + (Ce - C, 
C's) , as{mb)CF 



(Ca-Cs) [F^ + GK(srf) + Gi,(sb)] 



XG;,(./) + (C8 

f=u 



Nr. 



8G — CsG)GK,g 



+ C2G'K{sg) + [{G, - C,) + (Ce - C,)\ 5^G^(s/) + {Gsg - Csg)^^^^ 

f=u 



(32) 



9' 



where q takes u and c quarks, = (A'^^ — l)/(2Ai'c), and the loop functions Fx-, Gk, Gk. 
G'j^, and G'j^ g are given in Refs. [391 El- The internal quark mass in the penguin diagrams 
enters as s/ = m'j/ml^ The minus sign in front of (7^ (i = 3 — 6, 8G) comes from the parity 
of the final state as discussed in the previous section. 

By using above formulations, we estimate the SUSY contributions in the 6 — > d transition. 
In our calculations, we take /itan/3 to be 1 TeV and we set the MI parameters, |(5|'^)i3| = 
I ('^d'^)i3l ~ 10~^ from our previous works [3 [8]. We also assume that the magnitudes of the MI 
parameters {Sj^^)i3 and (5|^^)i3 are same but each phase is different. Thus, we parameterize 
the MI parameters as follows: 



|(5i«)i3|e--3 



\0d Jl3 - \{0d Ji3|e " . 



(33) 



''The A{B^ — )■ K^K'^) amplitude is explicitly presented in Refs. |39l|4T]. In our calculation, we neglect d 
(i = 8 — 10) since these Wilson coefficients are too small to contribute to the amplitude of — !• K^K^ in 
our model. 

'^The Cf (i = 3 - 6,8G) in Eq. ([32]) should be taken as the replacement Cf ^ [{VtbVt*^) / {VqbV*^)]C^ in 
Eq. Q. 
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Let us discuss the numerical analysis. In our calculations, we use the squark mass and the 
gluino mass as given in Eq. fl22|) . The present experimental data of BR(6 — )■ ^7) in Eq. fl25l) 
gives a constraint for the MI parameters as seen in Fig. [61 The SM contribution is larger 
than the SUSY one until |(5^^)i3| — 7 x 10~^, while the SUSY contribution dominates the 
6 — )■ 0^7 decay in the region of |((5|'^)i3| ^ 7 x 10~^. It is remarked that there is a lower bound 
of the branching ratio around 5 x 10^^. 

In Fig. [71 we show the allowed region of (5^^) 13 within 90% CL. of BR(6 — )■ ^7). It is 
found that any value of the phase is allowed in |(5^'^)i3| ^ 5 x 10~^. The upper bound of the 
MI parameter is at |(5|'^)i3| — 2 x 10"^ around the specific 9^^, 7i/2. 

By using this allowed region of (5^^)i3, we can predict the direct CP asymmetry A^f^"^^ 
and time dependent CP asymmetries Sp^, S^oj^o, and C^o^o. In Fig. [HI we show the predicted 
direct CP asymmetry A^^'''^ versus |(5^^)i3|- Here the value at 
Acp'^^(SM) ~ -0.09. Our prediction is -0.16 < Acp"^^ < 0.06 
future, we obtain an additional constraint of the MI parameters. 

In Fig. [H we show the prediction of Sp^ depending on |(5^^)i3|. The SM prediction is 



13 1 



is the SM one. 
If ^ct'^'*^ measured in the 



5p^(SM) ^ (2mrf/mfe) sin 0rf ~ 2.0 x 10" 



while the experimental data is 5'p^(exp) 



-0.8 ±0.7 [22]. In our prediction, the Sp^ reaches ±1 at 1(5^^)13! ^ 7 x 10"^ Therefore, the 



Spy is expected to be much larger than the SM prediction in the case of | (5 



LR\ I 



C(10- 



We expect the precise data to test our prediction in the future. 

In Figs. [TOl and [111 we show the predictions of the time dependent CP asymmetries 
Skoro and Cj^oro depending on |(5^'^)i3|, respectively. In the SM, one predicts 0.02 < 
S'^oi^o(SM) < 0.13 and —0.17 < C/^o^o(SM) < —0.15 [39J, while the experimental data are 
given as 5'^o^o(exp) = —0.8 ±0.5 and Cj^o^o(exp) = 0.0 ±0.4 [22], respectively. The present 
experimental bounds do not give any additional constraints to Fig. [71 However, more precise 
experimental data provide intensive constraints for MI parameters. 



3.0 

S 2.0 
? 1.5 

e 1.0 

m 0.5 
0.0 



0.0 0.5 1.0 1.5 2.0 2.5 

\i6]i%3\xlO^ 




|((5|/'^)i3|xl02 



Figure 6: The predicted region on the 
|(5|^)i3|-BR(6 ^ d-f) plane. The red 
solid and two red dotted lines denote the 
best fit value, upper and lower bounds of 
the experimental data with 90% C.L., re- 
spectively. 



Figure 7: The predicted region on the 
i^d^) 13^^13^ plane. The experimental con- 
straint of BR(6 — )■ d'-f) is taken account. 
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0.0 0.5 1.0 1.5 2.0 2.5 



Figure 8: The predicted direct CP asym- 
metry ^c'p'^'^ versus KS^j^^^ 
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Figure 10: The predicted time dependent 
CP asymmetry S^oxo versus |(5|'^)i3|. 
The red sohd and red dotted hues denote 
the best fit value and the experimental 
data with 90% C.L., respectively. 
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Figure 9: The predicted time dependent 
CP asymmetry S'p^ versus |(5^^)i3|. 




0.5 1.0 1.5 2.0 



2.5 



Figure 11: The predicted time dependent 
CP asymmetry C^o^o versus \{S^^)i3\. 
The red solid and two red dotted lines de- 
note the best fit value, upper and lower 
bounds of the experimental data with 
90% C.L., respectively. 



5 Summary 

We have discussed the contribution of the gluino-squark mediated flavor changing process to 
the CP violation in 6 — > s and h ^ d transitions taking account of recent experimental data. 
We have presented the allowed region of the MI parameters (5^^)23 and (5^'^) 13, which are 
constrained by the branching ratios of 6 — ?■ S7 and 5 — )■ ^7 decays. In addition, the time 
dependent CP asymmetries of — )■ J/tpKs, — )■ (pKs-, and B^ — )■ rj'K^ decays severely 
restrict the allowed region of the MI parameter, (5^'^) 23- These MI parameters [8^^)23 and 
(^d^)i3 s^ill allowed up to 2 x 10~^ for the squark and gluino masses of 1.5 TeV. If rriq ^ rrig 
increase, the bound of ((5^^)fc3 {k = 2, 1) is approximately rescaled as {S^^)k3 x mg/(1.5 TeV). 

By using these constraints, we predict the CP asymmetries of Bg — )■ (jxp, Bg — )■ ri'(j), and 
B^ — )■ K^K^ decays, as well as the CP asymmetries in 6 — )■ 57 and b ^ decays. We have 
summarized our results in Table [TJ It is remarked that the CP violation of the Bg — 00 
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Exp. 


SM 


our prediction 


BR(b s^) 


(3.53 ±0.24) X 10-^ 


(22J 




f3 15 ± 23) X 10"^ [371 


constraint 


BR(b d'y) 


(1.41 ±0.57) X 10"^ H 


11125] 


fl 54+^-^?) X 10"^ [251 


constraint 


^CP 


-0.008 ±0.029 [22 


] 




4 X 10-3 [sg 


-0.03 ~ 0.03 










-0.09 


-0.16 ~ 0.06 


c* 


0.679 ±0.020 [23] 






input 


constraint 


c 


0.741°;}^ [23] 






c 

- >^J/^pKs 


constraint 


c 


0.59 ±0.07 [23] 






c 

- '^J/V-A's 


constraint 


(t>s{Sj/^^ = sin 05 ) 


-0.004 ±0.166 ±0.054 [29] 


-U.Udbd_o.ooi5 [^ 


constraint 




-0.15 ±0.22 |22J 






= Sj/^^ 
= Sj/^^ 
0.04 [36] 


-0.2 ~ 0.4 
-0.5 ~ 0.4 
-0.4 ~ 0.2 


Spy 


-0.8 ± 0.7 






0.002 [36] 


-1 ~ 1 




-0.8 ±0.5 [22] 






0.02 ~ 0.13 [39] 


-1 ~ 1 




-0.0 ±0.4 [22] 






-0.17 0.15 [39] 


-1 ~ 1 



Table 1: Summary of the SM predictions, experimental values, and our predictions. 

decay is expected to be large owing to the squark flavor mixing. This prediction will be 
tested soon at LHCb. 
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